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In a chapter devoted to strain-energy methods Timoshenko 
and Goodter* have explained the principle of virtual work. They 
define virtual displacement in the ease of an elastic body as any 
small displacement compatible with the condition of continuity of 
the material and with the conditions for the displacements at the 
surface of the body, if auch conditions are prescribed. 

let us consider the case of plane stress distribution in 
@ plate. If from an orthotropic material a flat plate is cut 
parallel to a plane of elastic symmetry, it will have two perpendicu- 
lar axes of symmetry in the plane of the plate. Such a plate is 
Said to be orthotropic. 

Tf we denote by ~ and / the components of actual 
displacements due to the load and by /”~ and / the components 
of a virtual displacement from the loaded position of equilibriun, 
these latter components are arbitrary small quantities satisfying 
the conditions of continuity and the work done against the mutual 
actions between the particles is equal to the strain energy corres- 


ponding to the displacements. 
eae nian emer nemeneceeeneemen nme ee 
15, Tamoshenke and J, %, Goodier, Theory of 3 


Second Sdition (New Torente, and : 
Company, 1951), pig pt na Po dag i 
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The principle of virtual work states that the total work 
done during the virtual displacement of any system in eqailibrin 
mast vanish, Hence we have: 
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(2) 
where (/ means variation, X and ¥ are component body forces per 
unit volume of the plate, X and ¥ are components of body forces 
per unit area of the plate, and < and // are displacements in 
the x and y~directions, and where the first term in the bracket is 
the potential energy of deformation and the second and third terns 
together represent the potential energy of forces acting on the 
body if the potential energy of such forces for the unstressed 
condition 1s taken as sero, The complete expression in the 
brackets of equation (1) represents the total potential enerzy 
of the system, e can then say that the total potential energy 
of the system at the position of equilibriun is a maximum or a 
minimums and since, for stable equilibrium, it is always necessary 
to insist upon a positive work for any virtual displacement, the 
total potential energy of the system at this position is always 
& minima. 

For our problem we shall assume no body forces so that 

equation (1) becomes 
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For a plate of unit thickness the strain energy per unit 
volume Vo can be written as 
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In general the components of boundary forces per unit area 
of the plate can be expressed as 


XS My fmt, 
ES mty 4 Uy 
where 1 and m are the cosines of the angles between the normal 
and the axes x and y reapectively, 
For the problem which we wish to consider the boundary 
forces are to be exerted compressively normal to the boundaries 
y= ft caty cst we thewbing fees we capitis 
then X*o and ¥* ty (5) 
Yor orthotropic material the components of stress and strain 
are comected by the following relations s? 
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and for whieh /, and £, are Young's moduli in the x and y= 
directions, respectively, Poisson's ratio Cry ie the watto of 
the contraction parallel to the y-axis to the extension parallel 
to the x-axis associated with a tension parallel to the x-axis 
and similarly for fox oe emty “,, wane & 
rigidity associated with the directions of x and ys 

For an orthotropic plate like the one wo shall use, the 
NN fi Fs x, holds so that 2, = 2,, for our 


Since the equations of equilibrium are obtained independ- 
ently of the law connecting stresses and strains, they will be 
satisfied by a stress function 7 such that 
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With the knowledge of equations (3), (5)5 (6), (7) and 
(8) we can write equation (2) for our problem in the form 
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If we denote the expression in by // 5 we can write 
oV=o (20) 


Then a stress function taken in the form 
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can be determined from the minimum conditions 
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We shall apply this theory to the following problem, 
Ti. The Stress Distribution of a Wood Plate with 
Uniform Displacement 


Suppose we take a rectangular plate of Sitka spruce of 
unit thickness, length 2b, and width 2a, Referring to Figure ly 
let us prevent the corners A, B, C, D from any horizontal 
displacement, and place such a compressive load distribution 
(to be determined later) along sides AB and OD, only, as to cause 
the vertical displacements of all points along AB to be <k and 
along CD to be xk, where k is a positive constant, 

Stated forsally, the bouniary conditions which the 
stress function J must satisfy are: 
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at the corners (ab), (aynd)y (wagb)y (-a,~d), <= 0 
A suitable choles of the stress function ¢ for the 
problem under consideration is 
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whore Ay, (the k subscript on A, indicates the muber of 
constants Cg, that are included in each approximation) is an 
arbitrary constant to be adjusted to any desired total load on 
the sides y =» of the rectangle and /”,, is the Kronecker 
delta. 

Using equations (8) the stresses can readily be determined 
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We observe that /, > along /= +a 4n accord with the 
boundary conditions as givén by equations (13) 
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We also see that /,-° along /- +2 a9 should be true 
en ee ee (13). 
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Since 1, = oa » the first of equations (6) allows us 
to write 
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and, integrating with respect to 1, 
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To check our boundary conditions again we see thet << = 0 at 
the corners (a,b), (-a,d), faa (+ay=-b). 

ante! 5. 4 a » the second of equations 
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and, integrating with respect to y, we have 
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> a constant. 


Substituting from equations (15), (16), (17), (22) inte 


equation (10) gives 
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function 7 the greater the accuracy of our results, However, 
experience with past problems based upon the principle of virtual 
constants yields resulte of a high degree of accuracy, “We shall 
Limit our applied problem to the use of three constants in our 
discussions. 
ee 
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substitute inte equation (22), find a vatiaa as 
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Similarly, a8 @ third approximation, weing the three 
terns involving the coefficients C.5 Cz,» amd C.. 
in equation (13) we have 
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To dllustrate the problem for « particular case, suppose 
we take a plain-sawn plate of Sitka spruce with the grain of the 
wood parallel to the y-axis for which the elastic moduli are taken 
to be those given by Table 1, which follows, 


Oux 0.L6h 
fee 0.076 x 10° Ibs. per ea. ins 
Fy 1.679 x 10° Ibs. per sq. ins 


Ak 0.122 x 106 Ibs. per sas ins 
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Rie, hamise abl of Satins sprese se amusage wahaep 
Obtained from a large number of testa made at the U. 5, Forest 
Products Laboratory. 


Due to equations (7) the values given in Table I lead 
to the values given in Table Ii for By C2 a2, A 22) 


Ae: 


13.157 x 1076 ibs. per sq. ine 
0.59559 x 10° ibe. per eq. in. 

~ 0.27635 x 10°6 Ibe. por oqe ins 

~ 0.27635 x 1075 Ibs, per eqs in. 
8.9286 x 10° Ibs, per 69+ ins 


Using the values of Zmm given by Table II and r =k in 
equations (2h), (26), and (27) we arrive at the values of (47 
given in Table III, which follows, 
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TABLE III 
VALUES OF C,»,,,FOR SUCCESSIVE APPROXIMATIONS FOR SITKA 
SPRUCE 


Coy Cc 


First approximation # 2.2500 
Second approximation # 0.0366  #0.2hb70 
Third approximation #.0,03359 +f 0420620 = 0.03887 
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Sappose wa have a compression whose total load is P, 


where P is any positive constant. 
Gonsidering the case for the third approximation, we 


obtain the total load required by the formila 
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For our applied problem we must therefore tale 
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We are very much interested in the displacements «c 
and / along oe 
Pe 


letting {- +a = + / im equation (19) for the case 
where /c- 7 and substituting from equation (29) into 
equation (19) we have 
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ct ree cam bop 
(30) 


Ae. 


The results for displacement << are shown in Table IV 
and by @ solid Line in Figure 2a for the third approximation, 
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Por the third approximation we find that 
= 0.1592 4/0“ 
(33) 
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Hotties that. has the value 2 at two points along 
the boundary 1-2 whore y= 7b end ab wane -— 1b 
Using thres constants Cy; 5 C/2y Ca2 fram Table ITI, 
we have 
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The values of Uy eoog 4 = + / 09 chen oy ouation 
(37) are show in Table V and in Figure 3 by a solid line, 
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3/4 a - We 

5/8 a 7+ $660 

if2a — 4009 
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Ie a * 10,660 

f) ~ 12.620 
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“Yor rectangular plate of Sitka Spruce, 
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Tn the case of a square plate (a = b) equations (27) become 
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Solving equations (38) simultanecus ly for Case G2 
and, 15 we have 
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Again, to maintain a total load of Py we have 
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and substituting the values of 2,» C,. 5 C2~ given by equations 
(40) and a @b © 1, we find that 
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Now looking at the displacements again, we see that 
Mya ge FOE LIC Deft 1, 4004 Ce T#.00304] 
Or 
the values of < along ;= 7+. are tabulated in Table VI, 
which follows, and are plotted in Figure 4 by a solid line, 


TABIE VI 
VALUES OF DISPLACGMRAT zc AlOM X-*a (IN MULTTPI@S 
oF 10-6 p) * 
ha te 
b 0.0000 
7/8 tT 3140 
3/4» + Uel200 
5/8 + 21.7900 
ijf2b +t 20,0000 
3/8 b + 66460 
Yb 7 13,8700 
1/8 b 7 92.2600 
0 + 396100 
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The formila showing the nature of the load distribution 
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turns out to give 
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These values are tabulated in lable VII, which follows, 
and are plotted in Figure 5 by a solid line. 


eS FS ela TR 
x ic) = 
a f 34s 
78a 7 04638 
3/4 a 7 192080 
5/8 a 7 467920 
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Comparison of the resulte for the square plate with those 
for the rectangular plate shoms thet the results, in a general way, 
are the sam, 

ets quo nC aie ie; - rea 
taken after displacement under compression can be depended upon to be 
of the same general type regardless of the relative length and width 


of the plate. 
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For isotropic material 


Bes = es MD 2 
ie een 


(46) 


so that i 
ES <i LEE 

Aja ge Mae te ms P4337 ne 

Strectural steel is a good example of isotropic material. 
Suppose we consider a plate of such steel having the same dimensions 
and subject to the same bowslary conditions as we had for the wood 
plate and compare the displacement ~ along />ta for an 
isotropic material such as structural steel with the displacement 
A- @eng {- 12 for an orthotropic material such as Sitka 
spruce. 


Substitution from equations (46) into equations (27), using 
= 3 for structural stesly 2 = - + » gives 
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which, when solved simltanecusly, show that 
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To maintain a tetal load of compression P we tale 
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where the prime on the 7, signifies that we are dealing with 
the isotropic case, 
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whore, for structural steel, 8 #30 x 10° Ins, por square inch 
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Table VIII, which follows, shows the values of 
in miltiples of 107 p, 


ae = 


35 


TABIZ VITT 
% 
VALS OF ee), (Im woLtreizs oF 10°6 P) 
a er nn a ee nn ar enon F 
oe 
bv + 9000000 
7/8» + 9.002306 
3/4 b t Q,0077K4 
5/8 + 0012760 
if2b + 0.013880 
3/8 b + 0.009792 
Wav + 002252 
ife b = 0.004870 
r) — 0,007768 


The values of A. 4g ee Take VINE eo platted 
in Figure 2a as shown by a dotted Mine, Figure 2 shows in 
greatly exaggerated form the general shape of the plate during 
compression, 

They show that the general shape of the curve representing 
the dieplacement +. in the x-direction is the sam for both 
materials but that the displacements are very meh smaller for 
the isotropic material. 

The displacement / along 4. oem 

V mA | 


- — 2A “3 Cc me? — 
ie te tf yaa gi 2 (4 fate a Se /) (53) 


with substitutions from equations (46) and (48) giving 


tse pi Oradell Eg Raa = (54) 


This value is smaller than for the orthotropic material, 
as should be expected for the sam load P, 

The stress distribution along the top and bottem of the plate 
is given by 


/ 
yest ats 71 (/-C,)eo2h to (Cazes cot Yow 


sada 3 pe a 5 
a. LES SN ee 7] i 


Substituting the values of C,,,, from equations (48) 
and the values of 7,,,,, from equations (46) into equation (55) 
we have 
b).ee* (709 enh +: 4¥30en~ _, 2645+) P 
(56) 


Table IX, which follows, gives the values of {; along 
4th in miltiples of P, The valuse of Ky along 4 = + 
[SEE GREE EAB 
Lines 


aa 
7 907306 

8a > Os7912 
F 9073 
t 


i 


0.7075 


; 
[hoa ay 5 Eo 
bbeh 


or rectangular plate, isotropic material, 
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In the case of a square plate (a = b) of isotropic 
material using the three constants C2, 4 C). » Ci2 gives 


es Ca/- /6 00 Co t+ /€ 00 ae, Eat cr th 
—/6 00 Cay - SEFT Crt 6/72 Cur = “E70 (S7) 
SE 0O0Cr, FOI/FEfe ~ SSC Cie > EO 


which, when solved simultaneously, gives 


- #0,/7660O 


C2/ 
Rg al ne 


To maintain e totel coupression P we have 
o ame ‘ 
2, Wes we (4OntG, ¢Cr.4) A, = - 


As; = vate cy dc gk cs (59) 


Se Gaqtenet “< @long y=*+q is given by 
BA asl 
ee ok EZ 2 [Ceo Jen By a wig Go -)en Tp 


NK=La 


CI 4 El 2- Cay +4e2.)] ae) 


= £ (966AS CATE t, 07/300 EY 205) 178-60 
(62) 
The values of < along = +, ate tabulated in 


Table I, which follows, and plotted in Figure 4 as shown by a 
dotted line, 


TABLE X 
VALUES OF cc Alam X- 14 (iN MILI & 
106 p) * 


oar 

d 09,0000 
7/8 t 0.0197 
3/4 + 90,0518 
5f/@ + 0.0897 
172 b t 0.0763 
3/8 b + 0.0302 
if4b > 0.0418 
1fov ee 061066 
t 7 004325 


“Yor a square plate, isotropic material, 


The displacement in the yudirection is given by 


/ 
on WA, o rae ‘ss —/ 
a te Fe J2-T7 >/ sf ig Pe ) (62) 


ee Mae be 1/0“ 


ie (@) 


The values of hy ee 
equation 

(y PE Chik tiny ad +2,/387 pre 7 $002) P 
p98 


Similar to previous procedure. 
These values are tabulated in Table YI, which follows, 
and are plotted as shown in Figure 5 by a dotted line. 


TABIZ XI 
VALUES OF f ee 


a re 
7/8 a ZF 068234 
3/ha 7 1.9050 
5/8 @ 7 2aTT3O 
12a ZT 266390 
3/8 a 7. 12523 
1/4 a + 049048 
fa + 298480 

) + 346260 


Iv. Summary 

Before summarising we might mention the reason for 
setting up the boundary conditions in the manner set forth by 
equations (13), Dus te friction between machine heads and 
block, boundary conditions such as we have used might very easily 
be largely realized in actual practicon. 

Ag has been true of so many problems in the past, one might 
have expected the displacemnt < along (-— +a %o cause a simple 
bulging contour. instead, however, we get the unexpected contour 
shown in Figure 2b, 

Tt is very interesting to note that the average value of 
A g CMlled _X , does not necessarily occur at large distances 
from the ends (as one might expect), In fact, in some cases, it 
appears thet it occurs near the ends, This cautions against the 
common practice of placing the instruments for measuring Poisson's 
ratio near the middle, 

We also see that in the ease of isotropic material the 
unexpected contours still prevail but arc barely perceptible. 
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